Transport and the Order Parameter of Superconducting Sr 2 Ru04 
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Recent experiments make it appear more likely that the order parameter of the unconventional 
superconductor S^RuCU has a spin-triplet /-wave symmetry. We study ultrasonic absorption and 
thermal conductivity of superconducting S^RuCU and fit to the recent data for various /-wave 
candidates. It is shown that only / I 2_ B 2-wave symmetry can account qualitatively for the transport 
data. 



PACS numbers: 74.70.Pq, 74.20.Rp, 74.25.Ld 

The high-T c -analog compound Sr2Ru04 poses a fun- 
damental challenge to the study of unconventional super- 
conductivity Q . Clean samples become superconducting 
at about T c = 1.5K and the thermodynamics show that 
it is unconventional. The precise nature of its supercon- 
ducting order parameter remains elusive, however. The 
simplest spin triplet case is the time-reversal-symmetry- 
breaking odd-parity p-wave state which has an order pa- 
rameter d = z(k x +ik y ) P). Spin-triplet character is con- 
sistent with the Knight shift Q and neutron scattering Q 
measurements which indicate that the spin susceptibility 
shows no obvious suppression below T c . Moreover, /xSR 
experiment || confirms that the order parameter does in- 
deed break time-reversal symmetry. The central feature 
of this p-wave order parameter is that it exhibits finite 
gap all over the Fermi surface (FS). However, more re- 
cent thermodynamic and transport measurements at low 
temperatures have found powers laws in the temperature 
dependences. This includes specific heat Cy oc T 2 
penetration depth AA oc T 2 J7|, Ru spin-lattice relax- 
ation rate 1/T\ oc T 3 ||, ultrasonic attenuation a oc T 3 
H|0|, and thermal conductivity k oc T 2 These 
powers laws in T are consistent with line nodes in the 
order parameter. 

The combination of all experiments therefore favors 
the spin-triplet odd-parity /-wave symmetries which fea- 
ture both time-reversal symmetry breaking and line 
nodes. Maki and Yang jl3| have proposed an /- 
wave model with a superconducting order parameter 
d(T) = A(T) i k z (k x + iky) 2 (the so-called "3D" model) 
for Sr2Ru04. This 3D gap has line nodes where the 
basal plane (k z = 0) intersects the FS. It is similar to 
one of the leading candidates for the order parameter of 
UPt3 (known as the i?2« model in that hexagonal system 
fl4| ) . Graf and Balatsky [[l5| have proposed another can- 
didate /-wave order parameter d = A(T) z(k x -\-ik y )k x k y , 
which we will call the f xy model. This order parameter 
has line nodes perpendicular to the basal plane at where 
the planes k x = and k y = intersect the FS. This 
nodal structure is analogous to the p-wave model d = 
ATT) z(sin k x +i sin k y ) proposed by Miyake and Narikiyo 
fig . It is suggested that f xy -w&ve order parameter could 
arise from a pairing interaction mediated by antiferro- 
magnetic spin fluctuations with peak intensities at the 



incommensurate wavevectors q = (±2/37r, ±2/3tt). The 
latter has been confirmed by inelastic neutron scattering 
|fr7f . Finally there is the /^.^-wave order parameter 
d = A(T)z(k x + ik y )(k 2 — k 2 ) with line nodes where 
the planes k x — ±k y intersect the FS. This has been sug- 
gested [pi partly by analogy to high T c cuprates where a 
similar ~0r — y 2 (but with spin-singlet) gap structure has 
been conclusively established. It is worth noting that 
the above three models all account reasonably well for 
the specific heat and superfluid density data 18 19). 



The structure of all of these gaps in momentum space is 
rather complicated. Since all thermodynamic and trans- 
port quantities involve Fermi surface averages, it is not 
easy to see how they may be distinguished experimen- 
tally. In the study of heavy-fermion superconductors, it 
has been found that thermal conductivity and ultrasonic 
attenuation can be key tools in this situation. These 
quantities can be studied as a function of propagation 
direction and, in the case of the attenuation, also of po- 
larization dependence. The key point is that this di- 
rectional sensitivity can not only test for the presence of 
point or line nodes, but can also determine the positions 
of nodes and the orientations of nodal lines. These exper- 
iments [^(J offer the best chance to distinguish between 
order parameters that have different nodal structures but 
share similar thermodynamic properties. 

In this paper, we calculate the thermal conductivity 
and the transverse and longitudinal ultrasonic attenua- 
tion of Sr 2 RuC>4. We consider the three major candi- 
dates of /-wave symmetry mentioned above, namely the 
3D, f xy , and j x i_ y i models. It will be shown that only 
the f x 2_ y 2 model can account for the data even qualita- 
tively. The in-plane k 2 —k 2 symmetry offers an intriguing 
similarity to the high-T c order parameter. If this gap 
structure is indeed manifested in Sr2Ru04 it would sug- 
gest the presence of anti-ferromagnetic (AF) fluctuations 
in the ruthenate superconductors, perhaps peaked at or 
near q = (ir, tt) (such as one finds in the cuprates). The 
transport data thus suggests that AF fluctuation may 
also exist in ruthenates. 

Ultrasonic attenuation experiments in Sr2Ru04 have 
been done in the hydrodynamic regime where the sound 
frequency u> satisfies to << 1/r, where r is the electronic 
relaxation time fl(i|| . In this case the Boltzmann equa- 
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tion is appropriate, and the attenuation for sound 



propagation along the direction q 
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Here d denotes the dimensionality, F is the Fermi func- 
tion, and the angle brackets indicate an FS average. Note 
the second term inside the square bracket vanishes for 
the transverse case. Similarly the thermal conductivity 
in the zth direction when a temperature gradient is im- 
posed along the ith direction is computed using: 
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The FS of Sr2RuC>4 has been determined by de Haas- 
van Alphen experiments. It has three sheets denoted 
by a, (3, and 7 p2L These are in good agreement with 
LDA calculations p3fl. The bands are two-dimensional, 
a feature similar to cuprates. Among these three bands, 
the 7 band has the largest FS and has the highest den- 
sity of states (DOS) at the Fermi energy. For simplicity, 
we use a single cylindrical FS for f xy and / x 2_„2 models 
and a three-dimensional spherical FS for the 3D model. 
We shall comment on the effects of FS anisotropy below. 
We will assume that the momentum dependence of the 
gap is independent of temperature, which leaves only the 
overall magnitude A(T). The temperature dependence 
of the gap in Sr2Ru04 has not been directly measured. 
However, one may deduce the power law at low temper- 
atures from the pattern of nodes and near T c the power 
is given by mean field theory. For line nodes, a rea- 
sonable interpolation for the temperature dependence is 
A(T) = A(0)[1 - {T/Tcf] 1 ' 2 . This form has only one 
parameter A(0), reducing the arbitrariness in the fitting 
procedure. 

Inspection of Eqs. ([j]) and (||) shows that the only re- 
maining unknown is the relaxation time t s (T). Our 
procedure will be to use the thermal conductivity to de- 
termine the relaxation time and the magnitude of the 
ground state gap A(0), then use these to constrain the 
fits of the attenuation. 
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FIG. 1. Theoretical fits of thermal conductivity against the 
data of |Tl| in both the normal (upper) and superconducting 
(lower) states. The fits are equally well for the three candi- 
dates. The fitting parameters are given in the text. 

In Fig. [|, we first show the best fits against the ther- 
mal conductivity k x along the [100] direction in the nor- 
mal and superconducting states. The normal state is 
achieved experimentally by applying a sufficiently strong 
external magnetic field. The data are those of Tanatar et 
al. 0. In the normal state at 0.5T C <T <T C , the scat- 
tering rate can be well described by an empirical form: 
l/r„ = r (l + aT 2 /T 2 ) with T ~ 0.01ttT c and a = 0.15. 
That is, the inelastic scattering rate is observed to be 
roughly 15% of the elastic part at T c . For the effec- 
tive scattering rate in the superconducting state, we take 
the same temperature dependence as in the normal state. 
This is clearly appropriate for T <T C . At lower tempera- 
tures r«r c , due to the opening of the gap, the inelastic 
rate could decrease faster than T 2 . However, since the 
elastic scattering rate dominates over the inelastic part, 
the exact temperature dependence of inelastic rate is not 
so important. (This is one big difference between this 
system and the heavy fermion and high-T c materials.) 
Fitting to the superconducting thermal conductivity data 
then fixes the other parameter A(0). It turns out that all 
three theoretical candidates can fit the data for k x in the 
superconducting state. There is no noticeable difference 
in goodness-of-fit for the different models. The parame- 
ters determined by the fits are A(0) = 1.65T C for f x 2_ y 2 
model, A(0) = 3.3T C for f xy model, and A(0) = 4.05T C 
for 3D model. The difference in size between these var- 
ious values is not significant, since each gap form has a 
different relation between A(0) and the root-mean-square 
FS average gap. It is the latter quantity that is more 
important in determining the overall scattering strength. 

The parameters determined by the fit in Fig. [j] are 
now used to calculate the transverse and longitudinal ul- 
trasonic attenuation without any adjustable parameters 
(Fig. 0) for all three models. 
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FIG. 2. Theoretical calculations of ultrasonic attenuation 
for f x 2_ y 2, f xy , and 3D models (from top to bottom). The 
parameters used are identical to those in Fig. |l]. Experimental 
data is from Ref. JTo| |. 

As seen in Fig. the experiment exhibits the very 
clear feature that the transverse TlOO a is systematically 
higher than the longitudinal L100 a when scaled to the 
normal state values at T c . Note that the polarization 
vector for the transverse sound is also in the x — y plane. 
This sets a crucial constraint on theoretical models. 

Of the three candidate gaps, only the / x 2_ !) 2-wave has 
this property. The reason for this can be seen semi- 
quantitatively as follows. Inspection of Eq. ([!]) shows 
that one big difference between the transverse and longi- 
tudinal cases is a factor of k 2 k 2 : and (k 2 — l/d) 2 , respec- 
tively. The other factors restrict the integral to the FS, 
and, as a function of temperature, weight the points on 
the FS according to how small the gap is. At the lowest 
temperatures, only the nodal regions contribute. Thus 
we may think of the integral as being over an effective 
FS located near the nodes, with a polarization-dependent 
weighting factor. The 3D gap has an effective FS that 
is a circle at the equator k z = 0. For this geometri- 
cal object, we have {k 2 k 2 ) = ((k 2 - 1/d) 2 ) = 1/8 with 
d = 2. The angle brackets indicate an average over the 
effective FS. Thus the average for longitudinal polar- 
izations is about the same as that for transverse. For 
the f xy state, the effective FS is a set of two lines at 



k x = and k y = 0, and we have (k 2 k 2 ) = 0, while 

{(k 2 — l/cf) 2 ) is finite (d = 2), favoring longitudinal po- 
larization. In contrast, for the f x 2_ y 2 state, the effec- 
tive FS is at k x = ±k y . The averages therefore lead to 
(fc 2 fc 2 ) finite, while ((fc 2 - l/d) 2 ) = (d = 2 again), fa- 
voring transverse attenuation. These averages show that 
the j x i_ y i will naturally have more transverse attenua- 
tion than its competitors. Although we have plotted only 
one set of curves for each of the states, this is a general 
feature. The fact that the attenuation data favor f x 2_ y 2 
over other candidates was already noted by Lupien et al. 
]iot without calculation. 

Since the data in Figs. |l| and || are taken on differ- 
ent samples (T c = 1.44K and 1.37K, respectively) , one 
does not expect a quantitatively good fit in Fig. |2| us- 
ing the same parameters in Fig. [l]. At first glance, our 
results of a for f x 2_ y 2 and f xy waves are qualitatively 
similar to Graf and Balatsky's ]l5j when a resonant scat- 
tering limit (5 — 90°) is taken under the self-consistent 
T-matrix approximation. Lupien et al. p0| pointed out 
that while the sound attenuation data favors f x 2_ y 2 wave, 
the low T power law is much lower than the calculated 
result of Graf and Balatsky for the longitudinal case. In 
contrast, as shown in Fig. ^| our calculation based on a 
phenomenological input of scattering time gives satisfac- 
tory good fit for the longitudinal case including the low 
T power law. It also leads to better fit for the transverse 
case especially at T < T c , compared to Graf and Bal- 
atsky's almost linear in T result. The still rather large 
discrepancy between the theoretical fits and the exper- 
imental data (especially for the TlOO mode) could be 
due to the order parameter which is more involved than 
having the pure f x 2_ y 2 symmetry. This remains to be 
studied in more details. 

It is important to note one major difference between 
our phenomenological treatment and the self-consistent 
T-matrix approximation used by Graf and Balatsky. In 
the unitary limit, the usual T-matrix approximation for 
the impurity scattering will assume a strongly frequency 
dependent scattering rate 1 /T; mp ~ 1/w in the super- 
conducting state [ |15| . In contrast, we have adopted 
a constant elastic scattering rate (i.e., Tq). Our phe- 
nomenological input relies truly on the good fit to the 
thermal conductivity data in the superconducting state 
(see Fig. @). 

In-plane thermal conductivity measured by Izawa et 
al. Jl^] with a rotating in-plane magnetic field has shown 
a four- fold pattern with a minimum when H || [110]. 
This angle dependence is consistent with the j x 2_ y i and 
3D models but rules out the f xy model jl2). In terms 
of variation of the amplitude, Izawa et al. claimed that 
the state d = A(T)z(k x + ik v )(cos ck z + \b\), originally 
proposed by Hasegawa et al. p4[ , is more consistent with 
the data based on the approach in Refs. |ll|[l8| . Here c is 
the spacing between two adjacent Ru layers and |6| < 1. 
This model has a circular-like line node parallel to the 
basal plane with any k z satisfying —1 < k z c < 1. This 



3 



gap is symmetrical for any given parallel k plane and will 
lead to a sound attenuation very similar to that of the 3D 
model above. Thus it also does not produce the proper 
polarization dependence. 

Finally we comment on the effect of Fermi surface 
anisotropy. While in Sr2Ru04 the FS has multiple sheets 
and there seems to be a coupling between the a and (3 
bands, the analysis of transport data given here will not 
be significantly modified providing that the D±h symme- 
try of the crystal is preserved. Moreover, our result for 
the 3D model stands if the spherical FS is replaced by 
an ellipsoidal one that has circular projection parallel to 
the basal plane. As long as transport data is concerned, 
it is likely that interband effects are unimportant. 

In conclusion, the transport data of Sr2Ru04 presently 
give strong preference to the / a ,2„ y 2-wave order param- 
eters with perpendicular line nodes at k x = ±h y . A 
sharper comparison of theory and experiment will be pos- 
sible if the data for a and K is taken on the same sample. 
Furthermore, we suggest that the out-of-plane transverse 
sound attenuation measurement should be performed un- 
der the polarization q\x and e || z. This will shed more 
light on the information of the location of line nodes and 
the overall nodal structure. 
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